IRAS 16293: A "Magnetic" Tale of Two Cores 
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ABSTRACT 

We present polarization observations of the dust continuum emission from the young 
star forming region IRAS 16293. These observations of IRAS 16293, which is a binary 
system, were conducted by the Submillimeter Array (SMA) at an observing frequency 
of 341.5 GHz (A ~ 880/xm) and with high angular resolution (~2"-3"). We find that 
the large scale global direction of the field, which is perpendicular to the observed 
polarization, appears to be along the dust ridge where the emission peaks. On smaller 
scales we find that the field structure is significantly different for the two components 
of the binary. The first component, source A, shows a magnetic field structure which 
is "hourglass" shaped as predicted from theoretical models of low mass star formation 
in the presence of strong magnetic fields. However, the other component, source B, 
shows a relatively ordered magnetic field with no evidence of any deformation. We have 
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possibly detected a third younger outflow from source A as seen in the SiO emission 
which is in addition to the two well known powerful bipolar outflows in this kinematically 
active region. There is an observed decrease in polarization towards the center and this 
"polarization hole" is similar to decreases seen in other young star forming regions. Our 
calculations show that in IRAS 16293 the magnetic energy is stronger than the turbulent 
energy but is approximately similar to the centrifugal energy. There is considerable 
misalignment between the outflow direction and the magnetic field axis and this is 
roughly in agreement with model predictions where the magnetic energy is comparable 
to the centrifugal energy. In conjunction with other observations of the kinematics as 
determined from the outflow energetics and chemical differentiation we find that our 
results provide additional evidence to show that the two protostars appear to be in 
different stages during their evolution. 

Subject headings: ISM: individual (IRAS 16293-2422) - ISM: magnetic fields - polar- 
ization - stars: formation - submillimeter - techniques: polarimetric 



Introduction 



In the "classical picture" of star formation, magnetic fields are believed to strongly influence 
star formation activity in molecular clouds. They provide support to a clo ud against gravita- 



tional collapse a nd thus explain the low efficiency of the star formation process (|Mouschoviasll2001 
Shu et al.l 120071 ) . The process of ambipolar diffusion, in which magnetic flux is redistributed in 



the cloud, leads to the formation of a core that can no longer be magnetically supported and 
gravitational collapse sets in. In addition, the process of magnetic brakin g can help to remove 
angu lar momentum and slow down the rotation of the cloud as it collapses (jBasu fc: Mouschovias 
19941 ). In contrast, there have been a number of alternate theories which postulate that magnetic 
fields are relatively weak an d supersonic magnetohydrodynamic turbulence is the dominant process 
(|Mac Low &; Klessenl 12004! ) . Turbulence controls the evolutions of clouds, and cores form at the 
intersection of supersonic turbulent flows. Only a fraction of such cores become supercritica l and 
collapse begins to occur on a gravitational free-fall timescale (jBallesteros-Paredes et al.l 120071 ). 



Studying the morphology of the magnetic field in the interstellar medium (IS M) requires dif- 
ficult and sensitive observations of polarized radiation (jHildebrand Kirbyll2004l ). Spinning dust 
grains in the ISM become partially aligned with the magnetic field, generally with their long axes 
perpendicular to the field (jDavis &: Greensteinlll9 51). Nevertheless, the exact nature of the align- 
ment process is a matter of debate ( Lazarian 20071). but it appears that the do minant mechanism 
may actually be alignment by radiative torques (jHoang &; Lazarianll2008l . |2009| ). As a consequence 
of this alignment, the thermal dust emission is partially linearly polarized, with polarization direc- 
tion perpendicular to the magnetic field. The magnitude of the field can be determined indirectly by 
using the dispersion of the position angles of the field under the assumption of energy equipartition 
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between the kinetic and perturbed magnetic energies ( Chand rasekhar fe Fermilll953l ; lOstriker et al 
200ll : iHeitsch et aDhoOll : ICrutcher et aDhooi iFalceta-Goncalves et aliboosh . 



Observations of the large scale polarization distribution in molecular clouds from optical, infra- 
red, and submillimeter (submm) wavelengths usually s how an ordered patter n (IDotson et al.ll2000l : 
Pereyra fc Magalhaed 12004 : IPoidevin &: Bastienl 120061 : Matthews et al.1 12003 ) . T he analysis of po- 
larization data ha s showed that some dark clouds are magnetically dominated (jAlves et al.ll2008l : 
Hever et al.l 120081) while other giant molecular clouds are close to equipartition with turbulence 
([Novak et al.ll2009l ). However, o bservationally, the deb ate on whether turbulence or magnetic fields 
are dominant is still unresolved (jCrutcher et al.ll2009l ). Early and pioneering work with the BIMA 
millimeter (mm) array showed the dust (and molecular line) polarized emission can be traced 
at an angular resolution of few arcseconds, the scales where the core contr action process occurs 



(Rao et al 



1998 



Girart et al. 



1999 



Lai et al. 



2001 



20031 : ICortes et al.ll2005h . Currently, the Sub 



millimeter Array (SMA) is the only telescope that can detected the polarized emission at mm 



2009 



Tang et al 



2009 



30). 



In 



and submm wavelengths (jMarrone et al.1 120061 : iGirart et al.1 12006. 
the object NGC 1333 IRAS 4A (IRAS 4A here after), the SMA observations have revealed that 
the magnetic field configuration is consistent with theoretical models for the formation of solar 



type stars in which m agnetic fields play a much stronger role than turbulence (jGirart et al. 
Goncalves et al.ll2008h . 



2006 



The low-mass star forming region IRAS 16293 — 2422 (IRAS 16293 here after), located in 
the p Ophiuchi molecular clo ud complex, has been the focus of nume rous studies since the report 



of infall spectral signatures (jWalker et al 



1986 



Menten et al.1 119871 ) . Accurate observations us- 
i ng astrometric tech n iques show that the p Ophiuchi molecular cloud is at a distance ~ 120 pc 



(jKnude Hoe 



1998 



Loinard et al 



2008), while maser VLBI observations towards IRAS 16293 



suggest a distance of ~ 178 pc (jlmai et al.l 120071 ) . For this paper we have assumed a distance of 
d = 150 pc. IRAS 16293 is a Class protostellar system with a bolometric luminosity of 32 L , 
and is surrounded b y a compact, R env ~ 3000 AU, but relatively massive envelope, ~ 3.0 M , 
(jCorreia et al.ll2004l ). Higher resolution radio interferometric observations have revealed a double 
core separated by 5" (750 AU in the plane of the sky). The southern core is commonly referred 
to as source A and the northern one as source B. Both cores appear to have d i fferent physical 



and chemical properties (|Woottenlll989i : lEstalella et al.lll99ll : iBottinelli et al.ll2004l : IChandler et al 



20051 : iTakakuwa et al.ll2007l ). Despite its low luminosity, IRAS 16293 has a rich chemistry, with hot 



core l ike properties at scales of ~ 100 AU (IBlake et al 



200J; ICazaux et al.l 120031 ; iKuan et al.1 12004: iB isscho p et al 



lar outflow asso ciated with source A (| Walker et al.l Il988l : 



1994; 



Ceccarelli et al. 2000; Schoier et al 



2008 1) . There is a strong quadrupo- 



Mizuno et al.l Il990l ; IStark et al.l |200J; 



Yeh et al.l l2008l ) . The larger bipolar outflow is in the east- west direction (E-W outflow), with a 
position angle (PA) of 110°. The other outflow is in the no rtheast-southwest d irection (NE-SW 



outflow), with a PA o f 60° and shows copio us SiO emission (jHirano et al.ll200ll ). High resolution 



submm observations (jChandler et al.l 120051 ) have revealed multipl icity in source A an d which of 
these sources powering these two outflows is a matter of debate (ILoinard et al.l 120071 ) . Previous 



- 4 - 



measurements of the magnetic fiel d geometry in this ob j ect by observing polarized dust emission 
have yielded contradictory results ([Flett &: Murray! 1 1 99 ll ; iTamura et al.lll993l ; lAkeson &: Carlstrom 
19971 ). possibly because they were sensitivity limited. Thus, it is apparent that more sensitive and 
higher resolution observations of the magnetic field structure are needed. 

In § 2 we briefly describe the observations and the data reduction procedure. § 3 presents the 
results, § 4 the analysis, and § 5 discusses the possible scenarios for the observed magnetic field 
morphology in context with the information from observations already known from the literature. 



2. Observations and Data Analysis 



The observations were conducted in April 2006 (se e details i n Table!]) with the SMA0, which 
is located near the summit of Mauna Kea in Hawaii (|Ho et al.1 12004 ) . As can be seen from the 
Caltech Submillimeter Observatory (CSO) tau meter atmospheric opacities in Table[TJ the observing 
conditions on the first day were significantly better than those on the second day. Nevertheless, the 
atmospheric conditions were stable on both days resulting in minimal fluctuations in the antenna 
gains. The coordinates of the pointing center was at RA=16 h 32 m 22?9 and Dec.=-24°28'36'. / 00. 
The observing frequency was chosen to be located in the 345 GHz atmospheric window, which at 
the Mauna Kea site provides both optimal sensitivity and angular resolution. The SMA receivers 
operate in a double sideband mode with the two sidebands separated by ~10 GHz and the selected 
local oscillator frequency placed the lower and upper sideband central frequencies at 336.5 GHz and 
346.5 GHz respectively. The SMA correlator has a bandwidth of ~2 GHz which comprises of 24 
partially overlapping spectral windows. All the spectral windows had 128 channels, which provided 
a velocity resolution of 0.7 km s _1 . In addition to the continuum, the tuning frequency and the 
correlator were chosen to allow for simultaneous observation of the emission from the CO 3-2, 
SiO 8-7, and H 13 CO+ 4-3 spectral lines. 

Conducting polarimetric observations with interferometer arrays at mm and submm wave- 
lengths is challenging an d requires the use of s ome special techniques. A brief description of these 
techniques is provided in iMarrone et al.l (120061 ) and a more detailed d iscus sion of the meth o dolog y 



(both hardware and software aspects) is available in IMarron e (2006) and M arrone fe Rad ([2008) 
The data were reduced using the MIRIAD software package jwright fc Saultlll993h . The instru 



mental gains were calibrated by interspersing observations of IRAS 16293 (the target source) with 
observations of the quasars J1517 — 243 and J1622 — 297 which were used as gain calibrators. The 
instrumental spectral bandpass was calibrated from observations of the quasar 3C 273 and the 
absolute flux scale was determined from observations of Callisto. The single greatest factor that 
can corrupt the data is the intrinsic instrumental polarization which is commonly referred to as the 



1 The Submillimeter Array is a joint project between the Smithsonian Astrophysical Observatory and the Academia 
Sinica Institute of Astronomy and Astrophysics and is funded by the Smithsonian Institution and the Academia Sinica 
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"leakage". Therefore, the data were needed to be carefully calibrated to remove the effects due to 
this leakage. The primary task in MIRIAD which was used to solve for the leakage is GPCAL. For 
our observations we used 3C 273, whose intensity and polarization are strongly variable. In addition 
to the leakage, the task GPCAL also simultaneously solves for the polarization of the calibrator as 
well. At the epoch of our observation, we determined the intensity to be ~8.7 Jy and the linear 
polarization to be 1.0 ± 0.1% at a position angle of —34°. The leakages are different in each of the 
two sidebands as they differ in frequency by ~10 GHz. In the upper sideband, which is closer to 
the design frequency of the SMA polarimetry system, the measured leakages were approximately 
1%, while the lower sideband leakages were between 2 and 3%. These leakages were measured to 
an accuracy of 0.2% or better. 

The data from the source of interest, IRAS 16293, were then corrected for the leakages. The 
continuum emission in each of the two sidebands is contaminated by the emission from various 
spectral lines. The strongest spectral line emission that we detect is the one from the CO 3-2 
transition located at a rest frequency of 345.796 GHz. In addition, there is spectral line emission 
from various molecules such as SiO (8-7 transition) and H 13 CO + (4-3 transition), as well as SO2 
and CH3OH. There was also emission from other molecular transitions, but these were significantly 
weaker. The CO and SiO spectral lines are good tracers of the outflow activity in the earliest 
stages of star formation, whereas the H 13 CO + is a good tracer of the dense circumstellar gas. 
To create a pseudo-continuum channel for each sideband the contribution from the spectral lines 
was removed and the data were then averaged over all the spectral channels. These spectral 
line free and polarization calibrated data were used to produce maps of the I, Q, and U Stokes 
parameters. These maps were then independently deconvolved using the CLEAN algorithm. Since 
the continuum emission (Stokes I) from this source is quite strong, the source visibilities could be 
self-calibrated. The gain solutions from the self-calibration were applied to all the Stokes visibilities 
(continuum and line emission ) . The Q and U maps w ere combined to produce maps of the debiased 



linear polarization intensity (jLeahv &: Fernini Ill989i ). the fractional polarization, and the position 



angle. Continuum maps were made in each sideband separately for each of the observing dates and 
these were identical within the limits of the noise. Combined maps were then obtained using with 
the two sidebands and the two observing dates. Stokes I, Q and U maps of CO 3-2, SiO 8-7, and 
H 13 CO + 4-3 were also obtained. For the spectral lines, we present only the Stokes I line emission, 
since no significant polarized line emission was detected. Table [2] lists the basic parameters of the 
resulting maps, including the frequency of the spectral lines or continuum, the channel resolution 
(for line observations), the resulting synthesized beam, and the root mean square (rms) noise of 
the maps. 
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Results 



3.1. Continuum Emission 



The Stokes I map of the total continuum flux density is shown in Figure [TJ The contin- 
uum emission, which arises from t he warm dust, is quite strong and is resolved into the two well 
known components, sources A & B (|Woottenlll989l ). Subarcsecond resolution observations at A ~ 1 
mm show that source A itself is comprised of two components, Aa and Ab, separated by 0.6" 
([Chandler et al.ll2005l ) with Aa located southwest of Ab. In addition, subarcsecond VLA obser- 
vations show that the c entimeter (cm) wave length emission from Aa can be further split into two 
additional components (jLoinard et al.l 120071 ). The peak of the Stokes I continuum is ~ 4.1 and 
~ 3.9 Jy beam^ 1 for source A and B respectively, and the integrated flux from both sources is ap- 
proximately 11.5 Jy. The absolute flux scale at the SMA is only determined to an accuracy of 5% 
and thus all the flux densities calculated by us are uncertain by this factor. Single dish observations 
with the JCMT at 850 ^m (14V 6 be am) measure a flux of 23.6 db 1.3 Jy within a radius of 20" of 
the two sources (jCorreia et al.ll2004l ). The SMA measurements appear to detect about half of the 
total flux. This is due to the fact that the SMA antennas have a primary beam size (diameter) 
of 34" at the chosen observing frequency and the sensitivity to extended structure degrades as the 
distance from the pointing center increases. This missing undetected flux must therefore arise from 
emission on size scales larger than ~ 15" which is about six times the synthesized beam size. 

A Gaussian fit to source B shows that is barely resolved (see Table [3]) , with a radius less than 
100 AU. This is consistent with previous higher angular resolution observations w hich indicate that 



the dust emission comes f rom a optically thick disk with an outer radius of 26 AU (jRodriguez et al 



20051 : lLoinard et al.ll2007l ). The measured flux of sour ce B is in agreement with the value measured 
at at frequency of 305 GHz by I Chandler et al.l (|2005l ) with higher angular resolution, after taking 
into account its spectral index. The Gaussian fit to source A shows that this source is more extended, 
and is resolved with a deconvolved scale of 330 x 180 AU, elongated in th e north-south d i rectio n 
(see Table [3]). The flux measured is ~ 1 Jy higher than that expected from lChandler et al.l (|2005l ). 
probably because they filter some emission in their higher angular resolution maps. In addition 
to the contribution from sources A and B, there is some contribution from more extended dust, 
mostly from the southeastern part of source A and also north of source A (east of source B). The 
total flux measured with the two Gaussian fits of is ~ 10.6 Jy (Table [3|), whereas the total flux 
measured with the SMA is 11.5 Jy. 



3.2. Dust Polarization 

The linearly polarized component of the emission can be obtained from maps of Stokes Q and 
U. Typically, this is quite small and is only a few percent of the Stokes I emission. The maps for 
Stokes Q and U are plotted in in the top and bottom panels of Figure HJ The peak (absolute) 
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values of Stokes Q and U are ~ 7 times the noise level of ~ 4 mJy beam -1 . Note that in contrast 
to Stokes I which is a positive quantity, Q and U can be negative. We then obtained the maps 
of the (debiased) linearly polarized flux density (P), the polarization position angle (#), and the 
fractional polarization (p) which is expressed as a percentage. The maps of the errors in P, p, and 
9 are obtained as well. The map of the polarized intensity, fractional polarization, and position 
angle overlaid on a map of the total intensity is shown in the top panel of Figure [2j The fractional 
polarization and position angle are only computed at points where the debiased polarized flux 
density is greater than 8 mJy beam" 1 (~ 2a). Table H] contains a listing of the polarizations 
measured at various locations on the map. The errors in fractional polarization and position 
angle depend inversely on the Stokes I flux density and the polarized flux density respectively. 
Consequently, the errors in the fractional polarizations are smaller in regions where the continuum 
flux density is higher, while the position angle errors are smaller where the polarized flux density 
is larger. 

From the map of the polarized emission (top panel of Figure [2]) we can see that the polarization 
structures and morphologies are considerably different for the two sources A and B. The position 
angles of the linear polarization around source A appear to be approximately in a "centrosymmetric" 
pattern. Such a p attern can also be produced when the polarization is caused due to scattering 



i_pai 

3). 



(jSilber et al.ll2000l ). The scattering cross-section of the dust grains is quite significant in the optical 
and near-IR bands. However, it is inversely proportional to A 4 and thus the scattered polarized 
radiation decreases rapidly as the wavelength increases. Therefore, the contribution from scattering 
is likely to be negligible at the much longer submm observations. The polarization that is observed 
by us with the SMA must therefore arise largely from the continuum dust emission. 

The polarization fraction for the emission from source B is higher than that of source A, 
as seen in Figure [3j Furthermore, this plot shows that the decrease in fractional polarization at 
larger values of the total intensity is greater for source A. This decrease in fractional polarization 
towards the center of the cores where the peak emission is higher indicates that there may be a 
"polarization hole" effect. This decrease or depolariz ation is also seen in other sources as well 



(|Schleuningill998l : Matthews et al.ll200ll : lLai et al.ll2002l ). This effect could be due to at least three 



possible factors. Firstly, the higher density and temperature towards the center of the sources can 
lead to a misalignment due to the effects of a higher collisional rate and hence lower the degree of 
polarization. Secondly, the grain properties towards the center could be different from the grain 
properties in the outer envelope with the grains in the center being less able to align with the field 
than the grains in the envelope. The difference in the decreases in sources A and B could possibly 
be due to differences in the grain growth and properties in the two sources. One other possibility, 
which our observations seem to hint at, is that this decrease is due to the significan t disparity in the 



morph ology of the magnetic field. This central depolarization was also suggested bv lMatthews et al 



(|200ll ) from their observations of the OMC-3 filaments in Orion A. Towards the center of source 
A, the magnetic field directions change significantly over small scales leading to spatial variability 
in the position angles of the polarized emission. The resolution of our observations is not sufficient 
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to resolve the small scale structure. This results in lower polarization at the center due to the 
averaging of small scale structure. 

The polarized emission from this source has been the target of a number of different observa- 
tions conducted both with single dish telescopes and interferometer arrays at a range of wavelengths 
(Table [5]). As mentioned earlier in the introduction, most of the early measurements were incon- 
sistent with each other, which was perhaps due to limited sensitivity. A more recent polarization 
map from the JCMT using the Submillimeter Common Use Bolometer Array (SCU BA) polarimeter 



show s that at the peak of the dust continuum emission, no polarization is detected ([Matthews et al 



2009). We can compare our observations by convolving our images by the beam equal to the res- 
olution of the JCMT. Since almost all of the area over which the polarized dust emission detected 
by the SMA is within the JCMT beam, this is approximately equal to the integrated polarization 
over our map. When this is done we obtain extremely low polarizations around 0.2% (Table [5]), 
and is in good agreement with the SCUBA measurements. 

Under the assumption that the grains are not spherical and are rotating about their short axis 
which is aligned with the magnetic field, the field structure can be obtained by rotating the position 
angle of the observed polarization by 90° (bottom panel of Figure [2]) . There is considerable spatial 
structure in the deduced magnetic field directions, with a large scale twisted magnetic field in a 
direction coinciding with a curve joining sources A and B. In addition, source A shows that the 
lines are deformed with an "hourglass" like structure. This is the second such sensitive detection 
of this type of "hourglass" structure toward s a region of l ow ma ss star formation. The first was the 



young stellar object IRAS 4A reported by iGirart et al.1 ([200a). In contrast, source B shows very 
little variation in magnetic field structure. 



3.3. Molecular Lines 



3.3.1. H lz CO + 4-3 



Figure [5] shows the channel maps of the H 13 CO + 4-3 emission while Figure [5] shows the 
integrated emission and the velocity field. Both figures show that the H 13 CO + emission arises 
roughly extended in the north-south direction over approximately 17" x 10" (2500 x 1500 AU) and 
centered on source A. The strongest emission is offset by a few arcseconds to the north and south of 
source A. Indeed, the H 13 CO + integrated emission presents a relative minimum at the position of 
source A. Source B appears to be devoid of the H 13 CO + emission. As shown in Figures[Hand[5j there 
are clear signs of a velocity gradient along the major axis of the H 13 CO + structure (north-south 
direction) of about ~ 0.31 km s" 1 arcsec -1 , which translates to a physical scale of 430 km s _1 pc _1 , 
or an angular yelocity, 4.5 x 10~ 12 s - 1 . This is almost one order of magnitude higher that the 
value found by iNaravanan et al.l (|l998l ) from single-dish observations of the IRAS 16293 core. It 
is possible that this dis crepancy may be due to the different scales that are being probed. The 
Narayanan et al.l (jl998l ) scale-sizes are ~39" which is more than twice the maximum structures 
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being mapped by us. The north-south velocity gradient is also observed in the H2CO 4^3-3^2 line, 
although this line has a slightly different mo rphology: its emission arises from the two sources, A 
and B, and it is extended south of source A ([Chandler et al.l l2005). 



3.3.2. SiO 8-7 

Figure [6] shows the channel maps of the SiO 8-7 emission, which extends in velocity over 
10 km s~ 1 and is split into three main condensations. The most prominent one is observed just 
southeast of source A (SE condensation), and is detected in all the displayed channels, although it 
is brighter at the redshifted velocities. A Gaussian fit to the emission from this condensation in the 
^lsr = 6.94 km s _1 velocity channel gives a position angle of 145°, which is similar to the overall 
orientation of this condensation with respect to source A (see Figure [6]). The second condensation 
appears to arise northwest of source A (and southeast of source B) at systemic and blueshifted 
velocities (NW-1 condensation). At the ^lsr of 5.5 km s _1 the emission appears to break up into 
two condensations (also partially observed at the 4.14 km s _1 velocity channel). The peak intensity 
of the brightest channels appears to be located in the same axis as the one formed between SE 
condensation and source A. Indeed, the line that connects the SE and NW condensation passes 
closer to Ab than Aa. The third condensation is located about 6" northwest of source B (NW— 2 
condensation) . 



3.3.3. CO 3-2 

The high angular resolution maps of the CO emissio n (for the 2-1 an d 3-2 lines) from SMA 



observations have been already reported in the literature (|Yeh et al.ll2008l ). The emission is quite 
extended in the east- west direction (E-W outflow). In the visibility domain of our CO 3-2 data 
set, this shows up as a steep increase of the CO flux for visibilities with a radius, r UtV , shorter than 
20 kA. Figure [7] shows the channel maps of the CO 3-2 obtained by excluding the visibilities with 
r UyV < 20 kA. The most prominent emission comes from the E-W outflow. The blueshifted eastern 
lobe is more collimated than the redshifted western lobe, which has an open shell structure. The 
position of peak intensity of the brightest clumps (both eastern and western) appear to be well 
aligned, crossing source A, with a PA~ 105°. The SiO outflow is also traced by the CO 3-2 emission. 
The SiO SE condensation is detected in the lowest redshifted velocity channels (vlsr = 9.0 and 
11.1 km s _1 ). There is also blueshifted CO emission apparently associated with the SiO NW 1 
condensation, although the emission appears to be slightly displaced to the west, with the emission 
even being more bent to the west of source B. 
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4. Analysis 

4.1. Distribution of Polarization Position Angles 

The average values of the polarization PA around source A and B are 95° and 24° respectively. 
The distribution of the residual values of the polarization PA after subtraction of the average values 
around source A and B are shown in the bottom panel and the middle panel of Figure [8] respectively. 
This distribution shows that source B has a more uniform pattern than source A, which is likely 
due to the deformed morphology in the east-west direction (see Figure [2|) . In spite of the magnetic 
field structure, the analysis of the dispersion of the position angle as a function of displacement 
through the "structure function" can provide an indirect measurement of the turbulent to magnetic 



energy ratio (jHildebrand et al.ll2009l ). Since the range of scales where the polarization is detected is 



small, we pr efer to fit the magn etic field morphology around source A in a similar way as was done 



in IRAS 4 A (iGirart et al.ll2006l ). As a first approximation of the field geometry we fit the magnetic 
field vectors around source A only, excluding those around source B, with a set of parabolas. The 
results obtained with this method are in fair agreement with a more detailed analysis using specific 



theoretical magnetic field geometries (jGoncalves et al.ll2008l ). The parabolic functions used are of 
the type: 

(y - 2/0 ) = 9i + 9i C{x - x ) 2 (1) 

where yo and xq is the center of symmetry of the magnetic field configuration. The quantity gi is 
a parameter which depends on the curve that is selected and QiC represents the quadratic term. 
Note that the main magnetic field direction is along the x axis, and that for the observed magnetic 
field morphology, the y axis is going to be close to the right ascension axis. We used the \ 2 
method with yo, xq, C and the position angle of the main direction of the magnetic field, 6b as 
free parameters. The best fit solution obtained is with xq = — 0710, yo = —0^62, C = 0.25 and 
9b = 7.0°. Figure [9] shows the best solution for five sets of parabolas (top panel), and the values 
of the modeled magnetic vectors at the position of the observed magnetic field vectors (bottom 
panel). It is remarkable that this solution not only fits the magnetic field vectors around source A, 
but those around source B as well even though these were not used in the fitting! The distribution 
of the residuals, including the magnetic field vectors from source B, is shown in the top panel of 
Figure El The standard deviation of the residuals is 50 o b s = 13.9° ± 0.6°. The uncertainty of the 
polarization position angle is ag = 9.8° ± 3.8°, so the intrinsic dispersion is 56- in t = 9.8° ± 3.9°. 



4.2. Physical Parameters: Mass, Density and Column Density 

We can infer the total mass from the intensity of the dust emission using the following relation 
(Hildebrand 1983) 

Fd 2 
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where F is the flux, d is the distance to the source, k is the dust opacity, and B(T) is the Planck 
function. By fitting models to the single dish data obtained in the far infrared and submm regime, 
Correia et al.l (|2004l ) estimate the dust opacity law coefficient to be (3 ~ 1.8. From their results, we 
estimate the optical depth at 880 //m to be less than 0.18 at size scales greater than 150 AU (which 
is 95% of their fitted cloud size of 3000 AU). Thus, the emission is fairly opt ically thin over mos t 
of the cloud. This, however, is not the case towards the center of sou rce B ( Loinar d et al.ll2008l ). 



The dust temperature at the observed scales, is taken to be T~ 50 K (jCorreia et al.1 12004 ) . 



The integrated flux from sources A and B are given in Table [3j The calculated masses of 
so urces A and B are 0.33 M and 0.22 M Q respectively and these values are similar to ones derived 
bv lChandler et al.l (120051) . How ever, for source B, this mass is a factor of 2 lower than the value 
found by lRodrfguez et al.1 (j2005l ). and could possibly be due to the fact that the dust emission comes 
mainly from an optically thick disk. In order to determine the column density and the number 
density, we also need the total area used to estimate the flux density. From this area, we can define 
an equivalent radiu s of R = y/ 'Area /it. Following the same approach taken in the case of IRAS 4 A 
(jGirart et al.ll2006l ). we select the area of emission to be the entire region over which we can detect 
the continuum flux density (Stokes I). The equivalent radii, as calculated from the areas, are 4" and 
3" for sources A and B respectively. The column and volume densities source A and source B are 
given in Tableland is of the order 10 23 cm -2 and 10 7 cm -3 respectively. The combined total mass 
of this system according to our measurements is 0.55 M Q , with a column density of 6.4xl0 23 



cm 



and a volume density of 4.7xl0 7 cm" 3 . As discussed earlier, the total mass calculated by us is 
lower than the values determined from single dish measurements which are more sensitive to the 
emission on much larger scales. The mass of the envelope can be approximately determined if we 
assume that it contributes most of the missing flux of approximately 12 Jy. Using a value of 20 K 
for the envelope temperature, we obtain an envelope dust mass of approximately 1.9 M Q . This 
calculated mass is comparable to those derived from other observations of this source at a number 



(Walker et al. 


1990; 


Mezecer et al. 


1992; 



which is approximately 80% of the total mass of the system. This ratio is not ve ry different from 



that seen in other Class sources such as IRAS 4 A where its value is almost 90% (jJorgensen et al 



20071 ). The larger amount of envelope mass in IRAS 4A indicates that it is likely to be not as 
evolved as IRAS 16293. 

From the observed H 13 CO + velocity gradient, we can derive the dynamical mass needed for 
equilibrium between the gravitational and centrifugal forces: M^ yn = v^ ot R/G, where t> ro t is the 
rotation velocity and R is the radius of the flattened structure. For the measured values, R = 5'/2 
(780 AU) and v ro t = 3.31/sin 2 i km s _1 (i is the inclination angle of the rotation axis with respect to 
the line of sight), the dynamical mass is M& jn = 0.084 sin _2 i Mq. Assuming that the rotation axis 
and the outflow axis are parallel (as t hey appear to b e in projection), then the inclination angle 



derived from the outflow is i ~ 50-60° (|Yeh et al.ll2008l ). thereby M dyn ~ 0.11-0.14 M Q . Therefore 



the circumstellar mass around source A plus the mass already accreted onto the protostar is larger 
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than the dynamical mass, so this flattened structure is not stable, and likely is undergoing collapse, 
and has been inferred from spectral signatures of infall ([Chandler et al.l 120051 ; iRemiian Hollis 



20061 : iTakakuwa et alJl2007l ). 



4.3. Physical Parameters: Magnetic Field Properties 



The magnetic field strength can be estimat ed indirectly using two different methods: from th e 
modified Chandrasekhar- Fermi (C-F) equation ( Chandrasekhar Sz Fermi 1953 : Heitsch et al.|[200ll ): 
and from the curvature of the deformed, "hourglass" -like, field lines around source A. 



4-3.1. Modified Chandrasekhar- Fermi Method 



The modified C-F equation is: B = Q (5v/56) ^p, where 5v is the velocity dispersion along the 
line of sight, 56 is the intrinsic dispersion in the polarization position angles, p is the volume density, 
and Q is a dimensionless parameter that depends on the relative strengths of the magnetic field 
and the turbulence. We adopt a value of Q = 0.5, which is app ropriate for turbulent magnetized 
clouds with relatively strong fields, 56 < 25° (lOstriker et al.ll200ll ). We use the values derived in the 
previous section for the volume density and the intrinsic dispersion 56. The velocity dispersion is 
obtained from the H 13 CO + 4-3 spectral line emission as this spectral line approximately traces the 
same spatial scale as the polarization that is detected by us. From the intensity weighted velocity 
dispersion (second moment) map, the line of sight velocity dispersion is 5v ~ 0.35 km s _1 (i.e., a 
FWHM of 0.82 km s _1 ) in the regions where the emission is strong and is not affected by the strong 
velocity gradient se en in the north- south direction around source A. Then, using the modified C-F 
expression given in lLai et al.l ([2002]), we find that the component of the magnetic field strength on 



the plane of the sky is ~ 4.5 mG. 



4-3.2. Magnetic Field from Curvature of Field Lines 



The gravitational collapse of the cloud (neutral and ion particles) pulls the field lines into 
the canonical "hourglass" shape, producing a magnetic tension force resisting the collapse. This 
force, which is proportional to (B • A)B, can be approximately expressed as B 2 /R, where R is the 
radius of curvature. If the gravitational force is known, it is possible to estimate the magnetic field 
strength from the observ ed curvatur e of th e field lines using the following equation as derived from 
the expressions given by ISchleuninel (|1998l ) . 



B 



lmG 



R 



0.5pc 



D 



O.lpc 



M 



100 M(? 







n(H 2 



10 5 cm 3 



(3) 
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where D is the distance of a field line from the protostar. The radius of curvature of the magnetic 
field lines can be estimated from the family of fitted parabolic functions. For a parabola, y = a x 2 +b, 
the radius of curvature at the origin of the abscissa, is R = 1/(2 a). We selected the distance from the 
protostar to be 2" as this is approximately equal to our resolution along that direction. Furthermore, 
at larger distances along the center of source A, no polarization vectors are significantly detected. 
At this distance of 2" (300 AU), the radius of curvature of the field line around source A is ~ Of! 77 
(116 AU). The volume density measured previously is n(H2) — 4.9 x 10 7 cm -3 and the circumstellar 
mass for source A ~ 0.33 Mq. With these numbers, the magnetic field strength required for the 
observed curvature is about 3.5 mG, which is in reasonable agreement with the value estimated 
from the C-F method. 

Henceforth, we will estimate the relevant physical quantities using a field strength of ~4.5 mG. 
Using the estimated ayerage column density, we can calculate the mass-to-magnetic flux ratio 



(jMouschovias k, Spitzerl Il976l ) to be approximately equal to the critical value. This mass-to- 



magnetic flux ratio does not take into account the mass that is already accreted onto the protostars. 



From the modeling of the polarization pattern in IRAS 4A (jGoncalves et al.ll2008l ). the accreted 
mass in that source was similar to the mass as calculated from the dust. Therefore, our data sug- 
gests that the cores are approximately supercritical, which is in agreement with the fact that this is 
an active star forming site. The intrinsic dispersion of the polarization angles yields a turbulent to 
magnetic energy ratio of 0.32 ± 0.27, which suggests that the magnetic energy dominates over the 
turbulent energy. Finally, we can also estimate the relationship between the angular momentum 
and t he magnetic fie l ds, wh ich are also thought to play an important role in the dynamics of the col- 



lapse |Madiida_eLalJ (|2005l ). If the ratio between the angular velocity and the magnetic flux, uj/B, is 
larger than a critical value which can be expressed as, (uj / -E>) C rit = 3.21 x 10" 8 C' 1 yr' 1 fxG' 1 , then 
the angular momentum controls the collapse (c s is the sound speed in km s _1 ). If it is smaller, then 
the magnetic field dominates over the centrifugal forces. The measured magnetic field strength and 
the velocity gradient found in the H 13 CO + yields a ratio of oj/B = 2.8-4.8 x 10 -8 . This is slightly 
smaller than the critical value which is approximately 5.9 x 10 -8 (for eg = 0.54 km s _1 ). While the 
magnetic field appears to dominate energetically over the turbulence, it is comparable in magnitude 
to the centrifugal energy. 



5. Discussion 

5.1. Molecular Outflows, Kinematics, and the Evolutionary Stages of the 

IRAS 16293 Sources 



The IRAS 16293 region shows significant outflow activity with multiple outflows where the 
largest outflow is approximately in an east-west direction. O ur high-resolutio n CO 3-2 observations 
of this outflow are in agreement with previous observations (|Yeh et al.ll2008l ). that is, it is centered 
in source A, with an orientation of PA~ 105° and its morphology in the red lobe suggests that 
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this outflow is not very collimated (our maps were done excluding the short baselines thus probing 
extended structure) . The Spitzer near-IR images from the vibrationally excit ed H2 also trace sho ck 
structures within this outflow but at considerable distance from this object ([Padgett et al.ll2008l ). 



Our observations of the high-angular resolution C O 3-2 do not show any obvious traces of the 
extended NW-SE outflow, which was also not seen by lYeh et al.1 (120081 ). Therefore, it is possible 
that this is either a fossil outflow or that its powering source is presently in a very quiescent phase 
(with no strong accretion taking place). A third outflow has also been found in this region, as traced 
by the SiO 8-7 emission, centered in source A and with an orientation of PA~ 145°. This outflow 
is quite compact (less than 5" or 750 AU long) and very bright in the SiO line compared with its 
CO emission. The SiO 8-7 emi ssion in the outflow is expected to arise from dense (~ 10 7 cm -3 ) 



and hot (400 K) molecular gas (IHirano et al.1 120061 ) . All of this indicates that the SiO may trace a 
very young outflow, much younger than the two previously reported outflows. The powering source 
appears to certainly be within source A. The geometry of the outflow marginally suggests that 
source Ab may the powering source of this young outflow. Higher angular resolution observations 
are needed to attest the origin of this outflow. The structure of the NW blueshifted lobe at systemic 
velocities, which is split into two condensations with source B roughly in between, suggests that 
this outflow may be partially interacting with source B. 

The emission from this new SiO outflow can also be see n in the highest velocity channels of the 
combined SMA and JCMT maps of iTakakuwa et al.1 (120071 ) . They mention that the emission that 
they detect at high HCN 4-3 velocities is in agreement with the CO maps. The CO emission is quite 
complex because of contributions from both the EW and the new SE-NW outflow. The HCN 4-3 
emission at ambient velocities is possibly contaminated by this NW-SE outflow. Furthermore, the 
HCN 4-3 emission is also affected by the deep self-absorption and the high optical depths. Thus, 
it may not be a good tracer to try study the kinematics of the circumstellar /binary environment. 
In contrast, the emission from the SiO spectral line detected by us is much less contaminated and 
therefore, is a better tracer of this new outflow. 

Source B does not appear to show any active outflow and the only hint of any outflow activity 
may be the presence of a free-free component at very small scales (~ 15 AU) derived indirectly 
from the spectral index map at cm wavelengths (jLoinard et al.ll2007l ). This has been interpreted t o 



imply that source B has not yet started the phase of significant mass loss ([C handler et al 



Most of the emission comes from a very compact, optically thick disk ([Rodriguez et al.ll2005l ). with 



20051 ). 



a significantly smaller amount of material in the envelope around it when compared with other 
Class sources (see § 15 . 2|) . This suggests that most of the mass has already been accreted onto the 
disk and protostar. This raises the possibility that source B may not be a true Class protostar. 
The somewhat higher chemical rich ness of source A with respect to source B, may be due to the 
higher outflow activity of source A ([Chandler et al.ll2005l ). The lack of active accretion onto the 
disk from the circumstellar environment may be the cause of the non-existent outflow activity. It is 
possible that source B is (or was) the exciting source of the extended NW-SE outflow. Thus, source 
B may be a transition object between Class and Class I, or possibly even a Class I object (see 
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Ladal 119871 ; IWilking et al.l Il989l ; I Andre et al.l Il993l . for definitions and descriptions of the various 
classes of protostars). In this scenario, the uniform pattern of the magnetic field could be tracing 
the residual circumstellar envelope of source B. The narrower lin e widths associ ated with source 
B as compared to A have been interpreted as B being younger (jWoottenl 19891). If most of the 
material in source B comes from a mass ive disk in the plane of the sky ([Rodriguez et al.l 120051 : 



Chandler et al 



2005 



Loinard et al.1 120071 ). then kinematical motion in the disk will not broad the 



line width. 

Source A shows the more typical features of a Class protostar, with active and energetic 
molecular outflows surrounding a contracting core (as inferred from the "hourglass" morphology). 
This suggests that it is in an active accretion phase. The H 13 CO + 4-3 line traces dense molecular gas 
(its critical density is ~ 10 7 cm -1 ) and somewhat flattened structure in the north-south direction 
around source A. This flattened envelope is rotating with its axis in the east-west direction, which 
is in projection nearly para llel to the mo r e act ive and larger E-W. Interestingly, this rotational 
signature was first found by iMundy et al.l ()1990l ). who studied the kinematics in this region using 
the C 18 line. Their C 18 emission peaks roughly at source A and is extended with a PA=150°, 
and with a size that is not too different from the value determined from our H 13 CO + maps. 
Furthermore, they obtain the C 18 velocity gradient to be 2.4 km s _1 , whereas the value calculated 
from our H 13 CO + observations is ~3.0 km s _1 . This agreement suggests that both the C 18 and 
H 13 CO + spectral line emission seem to trace the same kinematical signatures of rotation! 

It is clear that the sources A and B are in different evolutionary stages and are likely to have 
different ages as well. Using su bmillimeter obsery ations of the source VLA 1623, which is also 
located in the Ophiuchus cloud, lAndre et al.l (|1993l ) find, based on calculations of the infall rate, 
that VLA 1623 is a very young Class source about 6000 years old. In contrast, recent surveys of 
star formation activity find that the Class sour ces in Ophiuchus a re at least three times older, with 
an estimated age of approximately 0.023 Myr ([Evans et al.l 12009). Even after accounting for this 
discrepancy, these Class protostars in Ophi uchus are ext r emely young when compared to similar 
sources in Perseus which are 0.32 Myr old. I Evans et al.1 (|2009l ) suggest that either the Class 
sources evolve rapidly into the Class I stage in Ophiuchus or that the transition from Class to 
Class I is not continuous. 

It is striking that the projected magnetic field configuration in the plane of sky is parallel to 
the elongation of the flattened structure and, thus, perpendicular to the rotation axis and the main 
outflow axis. One explanation for the lack of correlation between the outflow directions and the 
magnetic field could be due to the fact that the magnetic fields that we detect are mostly in the 
envelope while the disk ; from which the outfl ow originates, could be decoupled from the envelope 
and maybe precessing (jChandler et al.ll2005l ). A similar (but smaller) misalignment between the 
magnetic field direction and the outflow direction was also observed in IRAS 4A as well. Fur- 
thermore, single-dish observations of VLA 1623 als o find that the mag netic field is perpendicular 



to the outflow axis at scales of few thousand AU ([Holland et al.lll996l ). Most theoretical models 



that follow the collapse of magnetized and rotating cores assume that the magnetic field and the 
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rotation axis are aligned (jBaneriee fc Pudritzll2006l ; iGalli et al.ll2006l ). which, obviously, is not the 
case in IRAS 16293. Some recent simul ations have dealt with the situation of an oblique ma gnetic 
field with respect to the rotat ion axis (iM atsum oto Tomisakal 12004 ; iMachida et all 120061 ) . Ac- 
cording to the simulations of (jMatsumoto &; Tomisakal 12004 ). the angular momentum component 
perpendicular to the magnetic field axis is removed more rapidly than the parallel one, so the net 
effect is that the rotation axis becomes aligned wit h the magnetic field axis . This, however, does 
not a ppear to be the situation in our measurements. iMatsumoto et al.l (J2006J) have used the results 
from iMatsumoto &: Tomisakal (|2004l ) simulations in order to predict the polarization pattern of an 
initially oblique magnetic field (by 45°) in two scenarios, a strong magnetic field and a weak one. 
As stated before, magnetic braking at the scales relevant for the launch of the outflow (few tens of 
AU), induces an alignment of magnetic field and the rotation axis. Therefore, the outflow, which 
is already parallel to the field lines, maybe independent of the original field strength as well. On 
the other hand, the situation at larger scales (~ 1000 AU) is quite different. In the case where the 
magnetic energy density is strong compared to the centrifugal energy, the angular momentum is 
already aligned at these scales, thereby one should expect that the polarization observations shows 
the alignment of the B vectors with the outflow. If the centrifugal energy density is comparable 
to the magnetic energy density, the situation is significantly different, and the polarization obser- 
vations can show a projected magnetic field direction considerably misaligned with respect to the 
outflow. The extreme case occurs for a specific configuration (inclination angle with respect to the 
plane of sky direction), where the two axes can be almost perpendicular. It is noticeable that in 
this case, the polarization maps show a hint of an "hourglass" perpendicular to the outflow. 



5.2. Comparison of the Magnetic Field Structures in IRAS 16293 with IRAS 4A 



The magnetic field around source A shows the typical "hourglass" morphology that is expected 
from theoretical calculations. As discussed earlier, this is the second low mass star fo rming region 



wher e an "hourglass" configuration has been observed, the other one being IRAS 4A (jGirart et al 



2006). It is interesting to show that while there are some similarities, some significant disparities 
also exist as well: 



Both sources appear to be part of multiple star systems, with two dusty main components 
that have a sim i lar separation: 400 A U and 750 AU for IRAS 4A and IRAS 16293 respectively 
(IWoottenlll989l : lLooney et al.ll2000l ). The total bolometric luminosity, the total mass of the 
dense envelope, and the envelope radius surrounding the protostars are also not too different 



IRAS 16293 is only sligh tly more luminous but is less massive than IRAS 4A (jSandell et al 
199ll : ICorreia et al.ll2004l ) and powerful outflows originate from both objects. 



Another similarity is the contribution of the magnetic energy to the dynamics of the system . 
Since both sources appear to be contracting (IDi Francesco et al. ll200ll : IChandler et al.ll2005l ). 
their cores must be approximately supercritical. However the magnetic energy dominates over 
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the turbulent energy in both of the sources. 

• The biggest difference is that in IRAS 4A the "hourglass" morphology is detected in the 
massive circumbinary envelope, whereas in IRAS 16293, the "hourglass" is detected in the 
more compact circumstellar envelope around source A. In addition, towards IRAS 16293 most 
of the emission arises from the two circumstellar envelopes and from the circumstellar disks 
(particularly in source B), with only a small contribution from the circumbinary envelope, 
while in IRAS 4A the circumbinary contribution is significant. In order to quantitatively 
confirm this, we measured the average visibility amplitude in the 10-20 kA and 60-75 kA 
ranges for IRAS 16293, and in the 20-40 kA and 120-150 kA ranges for IRAS 4A (this source 
is about two times farther than IRAS 16293, so in order to trace the same physical scale 
the visibility range has to be twice as large). The amplitude ratio between the shortest and 
longest baselines are 2.1 ± 0.2 and 1.4 ± 0.1 for sources A and B, whereas it is significantly 
larger for IRAS 4A1 and IRAS 4A2, 4.1 ± 0.3 and 3.2 ± 0.2, respectively. 

• The total integrated fractional polarization from source A (~ 0.5%) is much lower than that 
obtained for IRAS 4A (~ 3.3%). One reason for this could be differences in grain properties 
that lead to low alignment efficiencies. Nevertheless, we cannot discard the possibility that 
the more compact emission detected in IRAS 16293 may also be responsible for this lower 
fractional polarization (beam smearing if there is complex unresolved magnetic structure). 

• A remarkable discrepancy is the difference in the orientations of the magnetic field axis and 
the outflow axis. For IRAS 4A both axes are not aligned, however the difference, ~ 40° is 
much less than that in IRAS 16293 where the main active outflow is nearly perpendicular to 
the main direction of the field (see § 15. ip . 



These differences suggest that IRAS 16293 is probably in a more evolved evolutionary stage 
and has more mass at smaller scales than IRAS 4A. The fraction of the total mass that is in the 
circumbinary envelope is larger in IRAS 4A. A significant amount of the mass has already fallen onto 
the circumstellar envelopes and the circumstellar disks of the two main protostellar components of 
the binary in IRAS 16293. The magnetic field configuration also suppor ts this scenario. It is also 
worth noting that the single-dish dust polarization maps in IRAS 4A (lAttard et al.ll2009l ) trace 
a uniform magneti c field which is in agreement with the higher angular r esolution SMA maps o f 



Girart et al.1 (|2006l ). whereas this is certainly not the case in IRAS 16293 (IMatthews et al.1 12009). 



Tentatively, this could be explained with IRAS 16293 being a more evolved region, where the 
outflow activity has put significant turbulent energy at the scales traced by the single-dish maps. 



5.3. The Magnetic Field Morphologies of IRAS 16293 on Various Sizescales 



IRAS 16293 is located close to the geometrical center of the L1689-northwest(NW) filament, 
which also contains the IRAS 16293-2422E core (which is approximately 2' to the east). This 
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filament extends about 10' (0.4 pc) in the NW-SE direction ([Nutter et al.l 120061 ). This is the 
direction that the overall dust emission appears to have from the SMA m aps and is also mor e 
apparent in the BIMA maps of the 2.7 mm dust emission at scales of 20" (jLoonev et al.ll2000l ). 
This large scale filament, which could conceivably be contained in a single magnetic flux tube, 
appears to bre ak up into at least tw o components towards IRAS 16293, sources A and B. Theoretical 
calculations by iMouschoviasI (ll99ll ) shows that the process of ambipolar diffusion can initiate single 
stage fragmentation along the length of a flux tube. The fragmentation, in the case of such a flux 
tube, is more likely to occur along its length which is parallel to the magnetic field, and is in 
agreement with our observations. Typically not more than three fragments are formed with each 
having a mass ~ 1M@. The fragmentation in magnetically subcritical clouds occurs when the 
hydromagnetic waves decay over length scales smaller than the Alfven length scale and this can 
occur when the densities are in the range of 10 3 -10 6 cm -3 . 

The SMA dust emission polarization observations conducted by us only sample the magnetic 
field on a relatively small scale (angular size < 15"). Comparable measurements with other in- 
struments (see § I3.2p were also similarly restricted in the sizescales probed. The magnetic field 
structure on much larger sca les were traced fr om observations of polarization due to dust absorp- 
tion at optical wavelengths (jVrba et al.lll976l ). They found that near the Ophiuchus cloud, the 
directions of polarization of the background starlight were approximately uniform. If so this must 
point to a strong and ordered magnetic field. Near-IR polarization observations, tha t probe deeper 



extinctions also show a direction that is similar to the optical polarization data (jWilking et al 



19791 : ISato et al.lll988l ). However, the nearest polarization vectors seen in absorption are at least 
> 1° away from IRAS 16293, which is roughly 3 pc projected in the plane of the sky and therefore 
the small and the large scale field structures cannot be easily connected. 

The ordered magnetic field structure observed by us indicates that the fields appear to be 
strong and are not significantly affected by turbulence on these small scales. However, this may not 
be the case on interme diate scales. The SCUBA polarimetry observations of the IRAS 16293 region 
(jMatthews et al.ll2009i ) containing sources A and B as well as source E show that the polarization 
position angles are considerably less ordered. This suggests that turbulence could possibly dominate 
on intermediate scales within the molecular cloud while magnetic fields regulate star formation 
activity on smaller scales where collapse signatures are seen in the mapped magnetic field topology. 
However, the sample of polarization vectors detected is small, and most of the detections have a 
low signal to noise (and consequently have greater position angle errors). Further higher sensitivity 
submm single-dish observations are needed in order to confirm if the dispersion observed is real. 

There have been some measurements of the strength of the line of sigh t component of the 
magn etic field (using the Zeeman effect in HI) in the p Ophiuchus cloud by iGoodman &: Heiles 
(| 19941 ). While they were able to measure the field strength significantly at a number of positions in 
this region, none of the detections are located in L1689. However, it is not clear whether this implies 
a low magnetic field strength in the low density gas component of the cloud, or a combination of 
projection (field close to the plane of the sky) and beam smearing (significant structure within the 
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beam gets averaged out) effects. 

6. Conclusions and Summary 

The installation of a polarimetry system on the SMA has enabled us to map the magnetic field 
structure in the ISM, especially in young star forming regions. Using this system, we have obtained 
high angular resolution and high sensitivity maps of the magnetic field structure in IRAS 16293 
through observations of the polarized dust continuum. These observations significantly improve 
on the earlier measurements which detected extremely low fractional polarization. Our detections 
indicate that those early attempts were limited in their sensitivity and angular resolution. At the 
same time, the net polarization over the entire region mapped by us is quite small (~ 0.2%) and this 
is indeed in agreement with past measurements. It is also apparent that the polarization fraction 
appears to decrease towards the center where the intensity of the source increases. However, it 
seems to taper off at sufficiently large values of the intensity especially in source B. This is possibly 
due to the effects of limited resolution towards the center where the peak of the emission occurs. 
Further higher resolution observations will be needed to determine if this relationship continues 
even at shorter distances from the center. 

The two sources A and B have significantly different magnetic field morphologies. In source 
A, our maps show that the magnetic field has the pinched "hourglass" shape that is expected from 
theoretical calculations. In contrast, the field lines in source B, appear to be quite uniform. Using 
the Chandrasekhar-Fermi method to obtain the field strength and the continuum dust emission to 
calculate the mass, we calculate the mass-to- flux ratio to be ~ 1. This implies that this object 
is in (or close to) the supercritical stage and the magnetic field is no longer able to prevent the 
collapse. In addition, the ordered field structures indicate that the magnetic energy likely exceeds 
the turbulent energy, but is comparable to the centrifugal energy. However, there are other reported 
measurements, that are not as sensitive, which appear to show that turbulence increases on more 
intermediate sizescales within the cloud. Nevertheless, the largest scales, which are probed in 
absorption polarimetry, once again show fairly ordered magnetic field structures. 

The evolutionary stages and the ages of the two sources are likely to be different and the SMA 
observations of the magnetic field structure also appear to confirm this view. The physical and 
chemical properties also differ in the two sources. In source B, most of the mass at scales of a 
few hundred AU has already been accreted onto a compact, optically thick disk, while source A 
still has a significant fraction of its mass in it's circumstellar envelope. Furthermore, the nature of 
the outflow activity is quite different in the two sources: source A shows significant activity and 
drives at least two powerful outflows and possibly a third compact outflow which we have detected 
in SiO emission, while source B shows no ongoing activity except for a possible fossil or remnant 
outflow. Furthermore, source A appears to have a richer chemical environment than source B. It 
is possible that source B may not be a true Class protostar and could possibly be a transitional 
object between Class and Class I. There appears to be a strong misalignment between the outflow 
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direction in source A and the magnetic field axis. This is in approximate agreement with theoretical 
model predictions when the magnetic energy is not significantly greater than the centrifugal energy. 

This is the second such low-mass star forming region in which it has been shown that magnetic 
fields appear to dominate over turbulence. However, the magnetic field morphology has been 
mapped in only a small number of such star forming regions. Improvements in telescope sensitivity 
to the polarized flux density will allow us to significantly expand this sample. With the increase in 
the number of objects studied, it will be possible to get a clearer and statistically significant picture 
and thereby address the all important question: Which process plays a dominant role in the star 
forming process — turbulence or magnetic fields? 
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Fig. 1. — Top panel: Color contour map of the Stokes U dust emission overlaid on the black contour 
map of the total dust emission (Stokes I). Stokes U blue (negative) and red (positive) contours are 
—5, —4, —3, —2, 2, 3, 4, 5, 6, and 7 times the rms noise of the map, 4 mJy beam -1 . Stokes I 
contours are 2, 7, 17, 37, 57, 77, and 97% of the peak in tensity (4.1 Jy beam -1 ). Crosses mark 
the position of Aa, Ab, and B sources (jLoinard et al.ll2007l ). The synthesized beam is shown in the 
bottom right. Bottom panel: Same as previous panels, but with Stokes Q in blue (negative) and 
red (positive) contours. 
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Fig. 2. — Top panel: Contour map of the total (Stokes I) dust emission overlaid on the grayscale 
image of the polarized dust intensity. The gray bars represent the polarization vectors. The 5% 
vector length is shown in the top left panel for comparison. Their length is proportional to the 
polarization fraction. Contours are 2, 4, 7, 11, 18, 28, 38, 48, 58, 68, 78, 88, and 98 percent of 
the peak intensity ( 4.1 Jy beam -1 ). Black crosses mark the position of Aa, Ab, and B sources 
(jLoinard et al.ll2007l ). The synthesized beam is shown in the bottom right. Bottom panel: Contour 
map of the total dust emission as in the top panel. The bars represent the magnetic field vectors. 
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Fig. 3. — The fractional polarization as a function of total intensity for each of the two sources A and 
B. The decrease in fractional polarization at large continuum Stokes I (also called a "p olarization 
hole") has also been seen at other polarization observations as well (|Schleuninglll998l ). For both 
sources the errorbars are also plotted and the fractional polarization errors are quite small compared 
to the value. 
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Fig. 4. — Contour channel maps of the H 13 CO + 4-3 line emission overlaid on the dust emission 
which is shown in the gray scale. The contour levels are —3, 3, 4, 5, 7 and 9 times the rms noise 
level, 0.42 Jy beam -1 . The i>lsr velocity of each channel is shown in the top left corner of the 
panels. 
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Fig. 5. — Top panel: Contour maps of the dust emission overlaid with the integrated intensity 
(zeroth moment) gray scale image of the H 13 CO + 4-3 line emission. Contours levels are the same 
as those in Fig. [2j Crosses mark the position of Aa, Ab, and B. Blue and red arrows show the 
directions of the E-W CO outflow and the NW-SE SiO outflow. Bottom panel: Contour maps of the 
integrated intensity (zeroth moment) overlaid with the intensity weighted velocity (first moment) 
color image of the H 13 CO + 4-3 line emission. Contour levels go from 5 to 95% the maximum value 
(8.7 Jy km s" 1 ) with steps of 15%. The units of the vertical bar is in km s^ 1 . Crosses mark the 
position of Aa, Ab, and B. 
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Fig. 6. — Channel contour maps of the SiO 8-7 line emission. For two of the channels, the dust 
continuum image in gray scale is overlaid. Contours are —3, 3, 4, 5, 7, 9, 11, 13 and 15 times the 
rms noise level, 0.35 Jy beam -1 . Crosses mark the position of Aa, Ab and B. The directions of 
the outflow is shown as red and blue arrows. The SE-NW arrows have a position angle of 145° and 
are centered on source Ab, and matches well with most of the SiO emission. The NE-SW arrows 
have a position angle of 105° and are centered on source Aa. 
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Fig. 7. — Channel maps of the CO 3-2 for the blueshifted (blue contours) and redshifted (red 
contours) emission. For the top-left panel the contours are 5, 10, 20, 30, . . . 120 times the rms 
noise level, 0.25 Jy beam -1 . For the other panels, the contours are 4, 6, 10, 15, 20, 30, 40, 50, 60, 
70 and 80 times the rms noise level. The fLSR velocity of each channel is indicated in the top left 
part of each panel. 
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Fig. 8. — From bottom to top: Position angle (PA) distribution of the polarization vectors around 
source A with respect to the average value, <PA>= 95.0°. PA distribution of the polarization 
vectors around source B with respect to the average value, <PA>= 24.0°. PA distribution of 
the residual from the parabolic fitting done to the source A data and applied to all the data (i.e. 
including source B), which is shown in Figure [9l 
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Fig. 9. — Top panel: Overlap of the 0.88 mm dust emission and the B-vectors (same as Figure [2]) 
with the best solution for the parabolic family of functions (blue solid lines). Bottom panel: Same 
as top panel, but the best family of functions are shown as B-vectors (in blue) at the same position 
as the measured values. 
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Table 1. Details of the SMA observations 



Date 


6 th April, 2006 


9 th April, 2009 


Opacity 51 


0.04 


0.1 


Array Configuration 


Compact 13 


Compact b 


Local Oscillator Frequency (GHz) 


341.5 


341.5 


Bandwidth per sideband (GHz) 


~2.0 


~2.0 


Spectral Channels per sideband 


3072 


3072 


Velocity resolution (km s -1 ) 


0.7 


0.7 



a Atmospheric attenuation at 225 GHz from the CSO tau meter 
b Provides projected baseline lengths from 7 meters to 70 meters 



Table 2. Continuum and spectral line parameters 







Synthesized Beam 


Spectral 


rms 




V 


HPBW 


PA 


Resolution 


Noise 


Observation 


(GHz) 


(arcsec) 


(deg) 


(km s -1 ) 


(mJy beam -1 ) 


Continuum 


341.5 


3.10 x 1.95 


-3.7 




4 a 


CO 3-2 


345.7960 


3.06 x 1.59 


-13.0 


2.11 


250 


SiO 8-7 


347.3307 


3.08 x 1.63 


-13.6 


1.40 


340 


H 13 CO+ 4-3 


346.9985 


3.16 x 1.84 


-7.9 


0.70 


420 



a Value for Stokes Q and U. The rms noise of the Stokes I map is affected by 
the limited SMA dynamic range, so it is higher, 14 mJy beam -1 
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Table 3. Properties of sources A and B 





Source A 


Source B 


R. A. (J2000) a 


16 /l 32 m 22!877 


16 /l 32 m 22!621 


Dec. (J2000) a 


-24°28'36'.'34 


-24°28'32744 


Deconvolved Size b 


2 / /29(2) x l'.'28(2) 


0"95(3) x / /59(6) 


Deconvolved P.A. b 


4° ± 1° 


85° ± 5° 


Peak Intensity (Jy beam" 


-l)a,c 4.io ±0.01 


3.89 ±0.01 


Flux Density (Jy) b,c 


6.14 ±0.04 


4.43 ± 0.03 


Flux Density (Jy) d,c 


6.85 ±0.05 


4.62 ± 0.04 


Equivalent Radius 


4" 


3" 


Fractional Pol. (%) c 


0.5±0.1 


1.0±0.1 


Pol. PA.(°) C 


104±6 


26±3 


Mass (M ) 


0.33 


0.22 


Volume Density (cm -3 ) 


4.9 x 10 7 


7.9 x 10 7 


Column Density (cm -2 ) 


5.9 x 10 23 


7.1 x 10 23 



a Estimated using Miriad's "MAXFIT" task. 

b Estimated from 2-dimensional Gaussian fittings to the image using AIPS's "IMFIT" task with 
a 0.5 Jy beam -1 cutoff to avoid the contribution from the weak extended component. Number in 
parenthesis give the uncertainty of the last decimal. 

c The absolute flux scale is accurate only upto 5%. 

d Estimated integrating the flux in the region used for polarization analysis, i.e. with a 
0.15 Jy beam -1 cutoff. 

c The fractional polarization and position angle were determined by integrating Stokes I, Q, and 
U over each source. 
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Table 4. Table of polarizations in IRAS 16293 



RA a 


Dec a 


I (Jy beam- 1 ) b 


P (mJy beam" 1 ) 


p(%) 


e (°) 


1.000 


-4.500 


0.342 


10 


3.0±1.3 


-51.1±12.5 


0.000 


-3.000 


1.218 


12 


1.0±0.4 


-75.0+11.0 


-1.000 


-3.000 


0.894 


14 


1.6±0.5 


-49.1±9.2 


0.000 


-1.500 


3.103 


16 


0.5±0.1 


-86.6+7.9 


-1.000 


-1.500 


2.402 


18 


0.8±0.2 


-63.0+7.2 


0.000 


0.000 


3.879 


17 


0.4±0.1 


-88.2+7.8 


-1.000 


0.000 


3.102 


19 


0.6±0.1 


86.4±6.8 


-2.000 


0.000 


0.981 


10 


1.1±0.5 


68.2±12.6 


2.000 


1.500 


0.192 


8 


4.2±2.4 


-40.6±16.0 


1.000 


1.500 


1.014 


18 


1.8±0.4 


-54.4±7.1 


0.000 


1.500 


2.180 


16 


0.7±0.2 


-84.0±8.0 


-1.000 


1.500 


1.893 


15 


0.8±0.2 


60.6±8.8 


-2.000 


1.500 


0.941 


19 


2.1±0.5 


48.4±6.7 


-3.000 


1.500 


1.069 


15 


1.4±0.4 


45.3±8.9 


-4.000 


1.500 


1.206 


15 


1.2±0.4 


23.3±8.8 


-5.000 


1.500 


0.413 


18 


4.4±1.1 


23.4±7.2 


1 nnn 

l.UUU 


^ nnn 


n ^7 


1 n 




-Oo.zinlo.O 


0.000 


3.000 


0.742 


9 


1.2±0.6 


-84.4±14.4 


-1.000 


3.000 


0.864 


8 


0.9±0.5 


56.0±16.5 


-2.000 


3.000 


1.089 


14 


1.3±0.4 


35.9±9.4 


-3.000 


3.000 


2.577 


29 


1.1±0.2 


30.5±4.5 


-4.000 


3.000 


3.459 


36 


1.0±0.1 


24.0±3.6 


-5.000 


3.000 


1.306 


25 


1.9±0.3 


14.1±5.3 


-3.000 


4.500 


2.089 


19 


0.9±0.2 


31.5±7.0 


-4.000 


4.500 


3.053 


28 


0.9±0.1 


30.9±4.6 


-5.000 


4.500 


1.259 


17 


1.4±0.4 


9.3±7.5 


-5.000 


6.000 


0.369 


11 


2.9±1.2 


-20.1±12.3 



a Offsets in arcseconds from the central pixel located at position with coordinates of 
RA=16 fe 32 m 22?9, and Dec=-24°28 / 36 / /00 

b Values of the Stokes I continuum flux density. The rms of the Stokes I continuum flux density 
map is ~17 mJy beam -1 
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c Values of the polarized flux density. The rms of the polarized flux density map is ~4 mJy 
beam -1 
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Table 5. Comparison of polarizations at different wavelengths 



Wavelength (mm) 


Beamsize (") 


p(%) 


on 


Reference 


2.84 


5.3x3.1 


0.7±0.5 


30±7 


Akeson et al. (1997) 


1.1 


19 


2.2±0.4 


135±5 


Tamura et al. (1995) 


0.88 


3.2x1.9 


0.17±0.01 


53±3 


This work a 


0.8 


11 


1.4±0.5 


62±11 


Flett & Murray (1991) 



a The polarization measurement is done by summing Stokes I, Q, and U in the entire map and 
then obtaining the polarized flux density, fractional polarization and position angles. 



